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ABSTRACT 

A new moving group comprising at least four Blue Horizontal Branch (BHB) stars 
is identified at {I, b) — (65°, 48°). The horizontal branch at go — 18.9 magnitude implies 
a distance of 50 kpc from the Sun. The heliocentric radial velocity is (Vr) = —157 ±4 
km s~^, corresponding to Vgsr = —10 km s~^; the dispersion in line-of-sight velocity 
is consistent with the instrumental errors for these stars. The mean metallicity of the 
moving group is [Fe/H]'^ —2.4, which is significantly more metal poor than the stellar 
spheroid. We estimate that the BHB stars in the outer halo have a mean metallicity 
of [Fe/H]=— 2.0, with a wide scatter and a distribution that does not change much as 
a function of distance from the Sun. We explore the systematics of SDSS DR7 surface 
gravity metallicity determinations for faint BHB stars, and present a technique for 
estimating the significance of clumps discovered in multidimensional data. This moving 
group cannot be distinguished in density, and highlights the need to collect many more 
spectra of Galactic stars to unravel the merger history of the Galaxy. 
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1 INTRODUCTION 

In recent years, stellar photometry primarily from the 
Sloan Digital Sky Survey (SDSS) has allowed astronomers 
to discover previously unknown spatial substructure in 
the Galactic spheroid. Other surveys such as 2MASS 
and QUEST have also been influential. Highlights in- 
clude the discove ry of new ti dal streams (Newberg et al.l 
20021: [Od cnkircti cn et al.l I2OO3I: fc rillmair & Johnson 200& 
Grillmair fc Diona tos 2006; Grill maidl2006l : lBelokurov et atj 
2006 : .Newberg et al.,.2009 | ): dwarf galaxies, s tar clusters, and 
transitio nal objects jWillrnan et al.l l2005a 3; Zucker et al.l 



2006a"^; 'Belokurov et al.' '20061, l2007l: llrwin et al.l bool 



Koposov ct al. 2009; Bclokuro Tet al.l 120091 '): and pre- 
viously unidentified spatial structure of unknown or 
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controversial idcntit 
2003.; .Ma icwski e t al 



jVivas et al.! I2OOII; lYannv et alj 



20041: iNewberg fc Yan ny 2005 



Belokurov et al.li2007lTjuric et al.ll2008l : lAn et al]l2009l ) 

The discovery of spatial substructure immediately sug- 
gests the need for spectroscopic follow-up to determine the 
orbital characteristics of tidal debris, the masses of bound 
structures, and the character of structures of unknown iden- 
tity. It has long been known that a spheroid that is formed 
through accretion will retain the record of its merger history 
much longer in the velocities of its component sta rs than it 
will in their spatial distribution (|Helmi et al.ll200^ ). For this 
reason early searches concentr ated on velocity to fin d "mov- 
ing groups" in the spheroid (|Maiewski et al.l 1 19961 ) rather 
than on the search for spatial substructure. 

Groups of spheroid stars that have similar spatial posi- 
tions and velocities (moving groups) are generally equated 
with tidal disruption of star clusters or dwarf galaxies, un- 
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like local disc "moving groups" that were once hy pothesized 
to be the r esult of disruption of star clusters ||Bo3I1934 
lEggenl ll958') but arc now thought to be the result of orbital 
resonances (Fam aev et al., 2008). 

Recently, searches for spheroid velocity substructure in 
SDSS data have concentrated on local (up t o 17.5 kpc from 
the Sun) metal-poor ma i n sequence stars ( Klement et al.l 
I2OO9I : ISmith et al.l l2009l : ISchlaufman et al.l boogf ). finding 
more than a doze n group s of stars with coherent velocities. 
ISchlaufman et all l|2009l ) estimate that there are 10^ cold 
substructures in the Milky Way's stellar halo. Blue horizon- 
tal branch (BHB) stars are particularly important for stud- 
ies of spheroid substruscture (e.g. Clewey & Kinman 2006) 
because they can be seen to large distances in the SDSS 
and because a large fraction of the SDSS stellar spectra are 
BHBs. In the future we need more complete spectroscopic 
and proper motion surveys designed to discover velocity sub- 
structure in the Milky Way's spheroid. 

In this paper we present evidence for a tidal mov- 
ing group of BHB stars, discovered in the SDSS and 
Sloan Extension for Galactic Understanding and Explo- 
ration (SEGUE; Yanny et al. 2009) spectroscopic survey. 
Additional evidence that the stars are part of a coherent 
structure comes from the unusually low metallicity of the 
stars in the moving group. We are unable to isolate this mov- 
ing group in density substructure, highlighting the power of 
velocity information to identify low density contrast sub- 
structure in the spheroid. 

Historically many of the spheroid substructures have 
been discovered by eye in incomplete datasets, rather than 
by a mechanical analysis of data with a well-understood 
background population. As the size of the substructures 
decreases, it becomes more critical to determine whether 
the observed structure could be a random fluctuation of the 
background. In this paper, we present a method (Appendix 
A) for estimating the probability that a random fluctuation 
could produce an observed "lump." 

As we work towards the fainter limits of the SDSS data, 
it becomes more important to understand how stellar pa- 
rameters derived from spectra are affected by a low signal- 
to-noise (S/N) ratio. The stars in the newly detected halo 
substructure have S/N ratio between 7 and 10, as measured 
from the ratio of fluctuations to continuum on the blue side 
of the spectrum. In the past we have successfully used spec- 
tra with S/N as low as 5 to measure radial velocities of F 
turnoff stars in the Virgo Overdensity (|Newberg et al.ir2007l ). 
In this paper, we show that although higher S/N is prefer- 
able, some information about metallicity and surface gravity 
can be gained for BHB stars with S/N ratios less than 10. 

Since the strongest evidence that these BHBs form a 
coherent group is that they have unusually low metallicity, 
even for the spheroid, we also explore the SEGUE metal- 
licity determinations for BHB stars in the outer spheroid. 
A great variety of previous authors have measured a large 
metallicity dispersion for spheroid stars, with average metal- 
licities somewhere in the —1.5 < [Fe/H] < —1.7 range 
l|Freemanll 19871 : iGilmore. Wvse. fc Kuiiken 1989). Those in- 
vestigations analysed glo bular cluster s ( Searlo fc Zinn 1978 ; 
IZinnl ll985|). RR Lvraes (ISahal Il985l : ISuntzeff et all ll99lF. 
K giants ( Morrison et al.l 2003 ) , and dwarfs ( Carney et al.l 
I1990I ). iNorrisI (|l986l ) found an average [Fe/H] of -1.67 for 
globular clusters in the outer halo, and an average metal- 



licity of —1.89 for fleld stars in the outer halo, again with 
a large metallicity dispersion that does not depend on dis- 
tance. The Besangon model of the Milky Way adopted an 
average [Fe/H] of — 1.7, with a = 0.25, for th e stell a r halo 
(|Robin et al.l I2OO0I ). Recently, ICarollo et all (|2007| . I2OI0I ) 
analysed the full space motions of > 10, 000 stars within 
4 kpc from the Sun, and found evidence that the ones that 
are kinematic members of the outer halo have an average 
metallicity of [Fe/H]=-2.2. 

Our analysis of the SDSS BHB stars shows that the 
metallicity of BHB stars in the Galactic spheroid does not 
change from go = 14.5 (6 kpc from the Sun) to go = 19.15 
(55 kpc from the Sun) at high Galactic latitudes. The mean 
measured metallicity of these stars is [Fe/H]=— 1.9, but anal- 
yses of globular clusters in the sample show that there could 
be a systematic shift in the BHB measured metallicities of 
a few tenths of a magnitude, and in fact [Fe/H]=— 2.0 is our 
best guess for the proper calibrated value. 



2 OBSERVATIONS AND DATA 

In order to search the stellar halo for co-moving 
groups of stars, we selected BHB stars from the sixth 
data release (DR6; Adelman-McCarthy et al. 2008) 
of the SDSS. This release contains both Legacy Sur- 
vey and SEGUE photometric data from 9583 square 
degrees of sky. More technical information for the 
SDSS 



survey can be found in 



York et al. 



Fukugita etahl ll 19961) : iGunn et all (Il998l): 



(2002 



. ISmi th et al.l 

IStoughton et al.l ||2002|) : lAbazaiian et al.l 




Pier et al.l (120031); Ivezic et al.r(|2004l ): iGunn et~aLl 
Tucker et all l|2006t ). 

We selected 8753 spectra with the colors of A stars 
from the SDSS DR6 Star database using the color cuts 
-0.3 < (ff - r)o < 0.0 and 0.8 < {u - g)o < 1.5, within 
the magnitude range 15 < (?o < 23, where the subscript 
indicates that the magnitudes h ave been correc t ed fo r red- 
dening and extinction using the ISchlegel et all (|l998l ) red- 
dening map. Within this color box, the stars that are bluer 
in (u—g)o tend to be high surface gravity blue straggler (BS) 
or main sequence stars, and those that are redder in (u — 3)0 
tend to be low surface gravity blue horizontal branch (BHB) 
stars. Using {u — g)o, {g — r)o color selection paramete rs as 
described in Figure 10 of lYannv. Newberg et al] (|2000| ). we 
divided the sample into 4630 candidate BHB stars and 4123 
candidate BS stars with spectra. 

We separated the stars in surface gravity using a pho- 
tometric technique, because in 2007 when we were selecting 
the stars, we were concerned that the S/N of the spectra 
would degrade more quickly than the S /N of the photome- 
try at the faint end of our sample. The S/N of the spec- 
troscopy increases over the whole magnitude range from 
15 < go < 20.5, while the photometric accuracy degrades 
only for u > 19. The photometric degradation affects halo 
BHB stars fainter than about go = 18. With more recent re- 
duction software (see below), we have found similar sensitiv- 
ity to surface gravity in photometric and spectroscopic indi- 
cators for faint SDSS spectra, but the spectroscopic separa- 
tion is much better than photometric separation for bright, 
high S/N spectra. The photometric selection technique does 
not introduce gradients in the selection effect with apparent 
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Figure 1. Identification of a moving group. The upper left and right plots show blue horizontal branch (BHB) stars with spectroscopic 
data in the region 50° < I < 80° and 40° < b < 55°, where the velocity is the line of sight velocity with respect to the Galactic standard 
of rest. The data in the two plots are the same, with the upper left plot marking two overdensities noticed on the upper right plot. 
Open circles in the lower left plot show the sky positions of the same BHB candidates with spectroscopic data in SDSS DR6. Filled 
circles correspond to stars that are in the co-moving stellar association, with the same stars plotted in the upper left and lower left 
plots. These stars have go magnitudes of 18.65 < go < 19.15 and line of sight velocities, with respect to the Galactic standard of rest, of 
—35 < Vgsr < 15 km s~^. The crossed circles show the positions of stars in a density at go = 16.5 found on the upper left plot. In the 
lower left figure, a higher fraction of the BHB stars have spectra near the globular cluster M13, at {I, b) = (59°, 41°). The box indicates 
the area of the newly discovered over-density. The smaller black dots in the lower left plot indicate the positions of photometrically 
selected BHB stars in Galactic coordinates. The lower right plot shows the positions of the photometric data from the lower left plot that 
has the go magnitude range of the moving group, 18.65 < go < 19.15. Note that the moving group is not detected from the photometric 
data alone. 



magnitude, and can be applied identically for stars with and 
without measured spectra. 

At bright magnitudes, the spectral S/N is high enough 
that the DR6 surface gravities are reliable. The SDSS spec- 
tra are R=1800, and the S/N of A star spectra varies from 
50 at go = 16 to 7 at go — 19. We tested the efficiency and 
completeness of the photometric selection technique using 
3996 SDSS spectra of stars with 16 < go < 17.5 and col- 
ors of A-type stars. Of the 1948 that had surface gravities 
consistent with BHB stars, 84% were classified as BHBs by 
the photometric technique. Of the 2048 spectra with high 
surface gravity (BS stars), 35% were classified as BHBs in 
the photometric technique. These numbers for completeness 
and contamination are similar to es timates that we made 
from examination of Figure 13 from lYannv. Newberg et al.l 
(|2000l ). and apply only to stars with u < 19. Fainter than 



this, the photometric errors in u will cause increased mixing 
of the populations. 



3 DETECTION OF A NEW MOVING GROUP 
IN THE GALACTIC HALO 

We searched our BHB catalog for stars that clustered 
in velocity, apparent magnitude, and Galactic coordi- 
nates. Our search was les s systematic than that used by 
IClewlev fc KinmanI (|2006l ). but allows us to find moving 
groups that are more extended, such as tidal streams. Be- 
cause they first searched for pairs of co-moving stars that 
were within 2 kpc of each other, they reduced their sensitiv- 
ity to tidal debris streams that could be spread over tens of 
kpc. 

We first divided the catalog into sky areas in {I, sin{b)) 
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that were 20° by 0.2 in the respective coordinates, and made 
plots of Vgsr vs. go magnitude. We used sin{b) so that each 
partition of the data would cover an equal area of the sky. 
Here, Vgsr refers to the line-of-sight velocity transformed 
to the Galactic standard of rest using the Solar motion o f 
ivx,VY,vz) = (10,225,7) kms"^ (|Dehnen fc Binnevll 19981 ). 
By examining the line-of-sight velocity as a function of mag- 
nitude, we found a number of co-moving groups of stars that 
were associated with known globular clusters, dwarf galax- 
ies, and tidal streams. Halo substructures that had been 
previously identified were not studied further. We concen- 
trated only on the most significant of the remaining BHB 
association; others were left for future studies. 

Figure 1 shows the Vgsr vs. go plot for the most sig- 
nificant co-moving group of BHB stars, along with the dis- 
tribution of these stars in {l,b). Although there is a second 
cluster of stars at go ~ 16.5 in the top left plot, we show in 
the lower left plot that these stars are dispersed in position 
on the sky so we did not explore that association further. 

The moving group we describe here has BHB stars at 
go — 18.9, agg = 0.1; and Vgsr = —10 km s~^, agsr = 10 km 
s~^. The dispersion is consistent with the instrumental cr- 
ror in the velocity of 19"^ magnitude BHB stars ( Sirko ct al. 
boo^ i. These stars have an average heliocentric radial veloc- 
ity of < (Vr) >= —157 ± 4 km s~^, which rules out any 
association of these stars with the disc population. In the 
direction {I, b) = (65°, 48°), the line-of-sight, heliocentric ve- 
locities expected for the thin disc, thick disc, and spheroid 
components are —2, —35, and —144 km s~^, respectively. 
The heliocentric radial velocity of our moving group is not 
far from that of the average spheroid star in that direction 
{Vgsr = for the spheroid and Vgsr = —10 km s~^ for the 
moving group), but the dispersion rules out an association 
with the general spheroid population. The lower right panel 
of Figure 1 shows that the association of horizontal branch 
stars could not be detected spatially; velocity information 
was required to identify the moving group. 

In the same area, at {l,b) ~ (59°, 46°), another 
cluster identified f r om S PSS BHB stars was found by 
IClewlev fc KinmanI (|2006l ). Although this cluster lies at 
roughly the same (l,b) angular position as the cluster iden- 
tified in this paper, the estimated distance to the Clewley et 
al. clump is 7.8 kpc, which is much closer to the Sun than 
our moving group, with a horizontal branch brighter than 
go = 16. 

The parameters of the seven BHB stars in the mov- 
ing group are presented in Table 1. The new moving group 
found in this study is located at {l,b) ~ (65°, 47.5°). The 
stars have apparent magnitude in the g filter of go ~ 18.9 
and velocities located within —35 < Vgsr < 15 km s~^, with 
an average velocity of < Vgsr >= —10 km s~^. A go ap- 
parent magnitude of 18.9 for BHB st ars, with an estimated 
absolute magnitude of Mg^ = 0.45 (|Newberg et al.l [2OO9I ') 
corresponds to a distance of ~ 50 kpc from the Sun. The 
magnitude dispersion of the seven stars is cr = 0.10. Note 
that the dispersion in apparent magnitude for these stars is 
agg = 0.1, and drops to agg = 0.03 if the two outliers with 
larger magnitudes (one of which has a high estimated sur- 
face gravity) are excluded. The magnitudes are consistent 
with a structure that has zero depth, but the depth could 
be as large as ct = 2 kpc. 



4 SURFACE GRAVITY ESTIMATES OF FAINT 
BHB STARS IN SDSS 

The SDSS DR7 SSPP pipeline (|Lee et al.ll2008al lb[) generates 
ten different measures of the surface gravity for each stellar 
spectrum, and additionally computes an "adopted" logg, 
which is determined from an analysis of the ten different 
methods. For stars with S/N less than 10, the adopted logg 
is set to an error code. 

We analyzed histograms of the distribution in logg of 
bright and faint A stars in our sample, looking for a method 
that showed two peaks: one at low surface gravity (BHB 
stars), and one at high surface gravity (BS or main sequence 
stars). Several of the methods appeared to separate bright 
A stars by surface gravity, and none appeared to separate 
the stars with S/N< 10 very well. We selected the logg9 
parameter, which is described in general terms in I Lee et al.l 
(|2008al '). in which it is referred to^as the "k24 grid" . The orig - 
inal reference to the method is iAUendc Prieto et al.l (|2006l ). 
where the method is described in great detail. The method 
finds the best match betwe en the SDSS s pectrum and a set 
of model atmospheres from lKurucj (| 19931 ). using the {g — r) 
color index and the normalized spectral fluxes in the wave- 
length range 4400 < A < 5500 A , using a resolving power 
of R — 1000. The surface gravities using this DR7 logg9 
estimator are tabulated in Table 1. Six of the seven candi- 
date members of the moving group have surface gravities 
less than 3.0, as expected for BHB stars. 

Subseq uent to ou r searc h for spheroid structure using 
BHB stars, IXue et all (|2008l ) published a list of SDSS DR6 
A stars along with measurements of Do. 2 a nd ,fm measure- 
ments. These are classical indices used by iPieij l| 19831 ) to 
classify BHB stars from the Hs line width and depth. Using 
this method for determining the surface gravities of A stars, 
only 2 of our candidate A stars appear to be low surface 
gravity. There are two effects that contribute to the small 
fraction of confirmed BHBs: (1) the performance of the D0.2 
indicator is degraded at low signal-to-noise, and (2) we dis- 
agree with the published Do. 2 values for two of the stars. 

In the top panels of Figure 2, we show the performance 
of the Dp 2 , .fm indicators for the set of 10, 224 stars from 
IXue et al] (|2008l ). and for the subset of 1176 of these stars 
for which the S/N of the spectrum is between 7 and 10. 
W hile the bri ght, high signal-to-noise spectra that were used 
bv lXue et al. I |2QQ8) are separable, the lower signal-to-noise 
spectra do not have two clear peaks in the diagram. We 
show histograms of Do. 2 for all of the stars between the two 
vertical lines (0.16 < fm < 0.24) in the top two panels of 
Figure 3. For all of the stars (most of which are high signal- 
to- noise), we see two peaks in the Do. 2 distribution. The low 
S/N sample shows no separation. 

We compare the separation using the properties of one 
line of the spectrum {HSD0.2 vs. fm) with the SSPP logg9 
surface gravities using photometry and the continuum levels 
of the spectra in the bottom panels of Figures 2 and 3. Figure 
2 shows that the logg9 indicator works about as well as the 
Do. 2 indicator for bright stars. Note that in the lower right 
panel, for stars with marginal S/N, the distribution of stars 
is much wider in the vertical direction, and the width in the 
horizontal axis is wider but not dramatically so. The lower 
panels of Figure 3 show histograms of logg9 for all of the 
stars between the two horizontal lines (15 < D0.2 < 40). 
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Table 1. Properties of BHBs in the Hercules Moving Group. 



ID 


RA 


DEC 


1 


b 


90 




(a - »■)() 




log(9) 


[Fe/H] 


S/N 


plate-mjd-fiber 


hh:mm:ss.s 


dd:mm:ss.s 


O 


O 


mag 


mag 


mag 


km s^i 


dex 


dex 




1336-52759-290 


16:08:33.7 


-1-38:49:19.7 


61.8 


47.6 


19.00 


1.22 


-0.21 


2.0 


2.75 


-2.99 


7.7 


1335-52824-1001 


16:06:59.7 


+40:20:39.5 


64.1 


47.8 


18.89 


1.17 


-0.16 


-4.4 


2.34 


-1.97 


8.5 


1335-52824- 167^ 


16:04:32.4 


+40:46:03.1 


64.8 


48.3 


18.85 


1.27 


-0.17 


-10.0 


2.91 


-2.86 


9.8 


1335-52824-362^ 


16:00:37.6 


+41:55:10.5 


66.6 


48.9 


18.82 


1.17 


-0.20 


-8.3 


2.83 


-2.63 


9.2 


1334-52764-540 


15:58:24.7 


+43:04:07.4 


68.4 


49.2 


19.10 


1.16 


-0.14 


-9.7 


3.39 


-2.10 


7.7 


1335-52824-008^ 


16:09:25.9 


+40:05:30.9 


63.7 


47.4 


18.86 


1.24 


-0.20 


-19.1 


2.92 


-2.27 


8.1 


0814-52443-046 


16:17:17.7 


+44:19:04.7 


69.7 


45.7 


18.89 


1.16 


-0.14 


-18.7 


2.62 


-2.42 


9.1 


1056-52764-637 


16:18:21.5 


+36:54:11.0 


59.1 


45.7 


18.90 


1.32 


-0.12 


-5.3 


2.21 


-2.99 


9.9 



1 Highly likely member of moving group. 



The brighter star data on the left shows two separate peaks 
for low and high surface gravity stars, while the marginal 
S/N data in the right panel does not show that separation. 
However, we can see from the correlation between D0.2 and 
logg9 in the lower right panel of Figure 2 that there is some 
information on the surface gravity of even these lower S/N 
spectra. 

Because we were surprised that the D0.2 numbers could 

be so high, even at low signal-to-noise, for stars that we 
expected were BHBs in a moving group, we independently 
computed new values for the width of the Hs at 80% of the 
continuum for the eight stars in Table 1. For six of them, we 
measured widths that were reasonably consistent with the 
published values. However, we measured a width of 23 for 
1335-52824-167 instead of 40, and a width of 23 for 1335- 
52824-008 instead of 33. With these corrections, six of the 
eight stars in Table 1 have Do. 2 widths less than 30, includ- 
ing all four of the stars that we will later show are highly 
likely members of the moving group. 



5 STATISTICAL SIGNIFICANCE OF THE 
MOVING GROUP 

Given that we were searching through the data by eye and 
looking for clusters, it is important to verify that the clus- 
tering of stars is not a chance coincidence. In this section 
we will estimate an upper bound for the probability that 
this clump could be a random fluctuation in the number of 
stars in a a small region of angle in the sky, magnitude, and 
line-of-sight velocity. 

Most of the sky that we searched for clumps was at 
high Galactic latitude; at low Galactic latitude the stellar 
population is different, the density of stars is higher, and 
the sky coverage is less uniform. So for statistical purposes 
we limited our sample to 6 > 35°. There are 2893 stars in 
our sample with b > 35°. The seven stars were discovered 
in a 19 square degree area of the sky (compared with 8800 
square degrees above b > 35°), with the requirement that 
this sky area be rectangular in Galactic latitude and longi- 
tude (61.8 < I < 69.7, 45.7 < 6 < 49.2). In comparison with 
the entire high latitude sample of BHB candidates, the seven 
stars in the clump span about 1/46*'' of the entire longitude 
range of the input data, and I/IO"' of the latitude range 
of the data in an equaJ-area plot. Prom Figure 1, we find 
that the range of velocities in the clump includes 15% of the 



sample stars, and the range of magnitudes in the clump is 
the same as 10% of those in the high latitude sample. 

Appendix A presents a statistical method to deter- 
mine the significance of a "clump" of data points in a d- 
dimensional search space. In this case we have four dimen- 
sions: Galactic longitude. Galactic latitude, radial velocity, 
and apparent magnitude. The statistical method measures 
the probability that random data will produce a clump of 
seven or more points within a four-dimensional box of the 
measured size. By using percentiles, we account for the non- 
uniformity of the data in the radial velocity and apparent 
magnitude dimensions. We also adjust the statistics, as ex- 
plained in Appendix A. 4. 2, for the periodic boundary condi- 
tion in Galactic longitude. What is important for the calcu- 
lation is that each of the dimensions be independent (sepa- 
rable, in the terms of Appendix A). In our case, the variables 
are not quite independent, especially Galactic coordinates. 
Because the density of A stars with SDSS spectra varies as 
a function of position in the sky, the relative percentiles in 
(/, 6) would be different if a clump of the same areal extent 
were discovered in a different part of the sky. We estimate 
that the local density in the region of the sky that the clump 
was found is 3.5 times higher than average. Therefore, fol- 
lowing the prescription of Appendix A4.4, we multiplied the 
longitude width of the clump by a factor of 3.5. 

It was not possible to account for every known peculiar- 
ity of the data. For instance, when we searched for clumps 
we divided the sky up into fixed (Z, b) boxes, so there are 
regions around the edge of each {I, b) box that were not as 
thoroughly searched. However, it is possible that if nothing 
was found with these limits we would have adjusted the box 
sizes and positions, so this oversight (which would underes- 
timate the significance of a clump) is probably justified. 

Given the percentiles of the 4-D search space listed 
above, we expect 0.32 stars in a typical box with these di- 
mensions. Using the algorithm in Appendix A, the number of 
clumps we should expect to find in this data that have seven 
or more stars in a 4-D box of these dimensions or smaller is 
0.96. Since we found one clump, this is very close to random. 
If we eliminate one star from the sample that is an outlier 
in {l,b), so that the area of sky is now 61.8° < I < 68.4°, 
47.4 < 6 < 49.2°, then we calculate that the E- value for 
finding six or more stars in a 4-D box with those (smaller) 
dimensions is 0.18, so the probability of finding this moving 
group is not more than one in six. Using only these dimen- 
sions, the clump of stars is not a significant detection. 
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logg9 logg9 

Figure 2. Lumin osity separation o f A colored stars into BHB and BS populations. Upper left panel reproduces the sample of 10,224 A 
colored stars from lXue et alj | |200S| ). showing that BHBs have lower H5 Z)o.2 widths than BSs. The upper right panel shows a low S/N 
subset of 1,176 stars (7 < S/N < 10). The candidate stream members from Table 1, which all have S/N in this range, are indicated with 
(red) crosses. The single Do. 2 method appears to be a less effective discriminant at S/N < 10. However, our independent measurement 
of Do, 2 for the moving group candidates resulted in significantly lower Hg widths for two of the eight stars (indicated be the vertical line 
below two of the crosses), so the effectiveness of the Do. 2 may depend on how it' s measurement is implemented for low signal-to-noise. 
The lower left panel shows D0.2 versus the 9**^ SSPP method (logg9)for the same lXue et al. 1 l l200Sf) sample. Note that the logg9 surface 
gravity discriminant achieves similar results to the D0.2 measure for bright stars. On the lower right, the lower S/N sample is shown. 
A cut at logg9 < 3.15 maintains the D0.2 separation into BHBs and BSs which appears effective even at S/N < 10. In this paper, we 
use an even more conservative cut at logg9< 3.0. Using this conservative selection criterion for BHB stars, 7/8 stream candidates (red 
crosses) are consistent with low surface gravity BHB stars. 



In order to claim a significant detection of the mov- 
ing group represented by these seven stars, it is necessary 
to compare their metallicities with the metallicities of the 
background population. Only 30% of the stars in the sam- 
ple have metallicities as low as those in the clump of seven 
stars. If we extend the statistical calculation to five dimen- 
sions, including metallicity, then the probability of finding 
seven stars in a region of this size and shape (or smaller), 
anywhere in the 5-D parameter space is not more than 1 
in 244 (E- value of 0.004), and the probability of finding six 
stars in a smaller angular region in the sky (in 5-D) is not 
more than 1 in 490 (E-value=0.002). These are highly sig- 
nificant detections. 



Since metallicity is critical to the detection of this sub- 
structure, we examine the SDSS DR7 metallicity measure- 
ments of these stars, and BHB stars in general, in the next 
two sections. In section 8, we use the information we have 
learned in the analysis of this moving group to re-select sim- 
ilar data from SDSS DR7, and make an even stronger argu- 
ment for the statistical significance of our result. 
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for the subset of the lXue et al.l (|2008|) sample with an HS flux minimum (as fraction of continum) of 0.16 < fm < 0.24. Upper right: We 
show the same histogram for the lower (7 < S/N < 10) signal-to- noise subsample. Stream candidates from Table 1 are indicated with 
a heavy curve. Note the two values at 40 A and 33 A should be moved lower (23 A ) based on inspection of individual spectra. Lower 
Left: We show the distribution of logg9 for most of the lXue et al.l ([2008') sample. Lower Right: We show the lower S/N subsample. Note 
the positions of the stream candidates (dark line), which are consistent with the low surface gravity BLIB population in 7/8 cases. Note 
the clear bimodality between BHBs (low Z?o.2 or low logg9) and BSs (high Z?o.2 or high logg9) for the full samples, and the absence of 
bimodality for the lower S/N samples. 



6 METALLICITIES OF THE MOVING GROUP 
STARS 

All seven of the candidate BHB stars in this moving group 
had unmeasured metallicities in SDSS DR6, but metallici- 
ties were assigned to these same spect ra in the much im- 
proved SSPP pipehne of SDSS DR7 (|Lee et all [ioOSal lbl) . 
Our metallicities come from the DR7, which became pub- 
lic in October 2008. The SDSS DR7 tabulates many mea- 
sures of the metallicities of stars, of which the most com- 
monly used is FEHA, the "adopted" metallicity, which is 
derived from a comparison of all methods used to measure 
metallicity. In this pape r we u se instead the metallicity of 
IWilhelm, Beers, fc Gravl l| 19991 ). hereafter WBG, which is 
specifically designed to measure the metallicities of BHB 
stars. The WBG metallicities of the seven BHB stars in the 
moving group are given in Table 1 and are shown in Figure 
4. The mean of the distribution is [Fe/H]= -2.46±0.14, and 
the width is cr = 0.4. 

We now show that the metallicities of the moving group 



stars are not consistent with being drawn at random from 
the stellar halo. We selected all stars in DR7 with the same 
photometric constraints as the stars in the moving group, 
and with b > 45° so that the BHB stars are likely to 
come from the same component (the spheroid), and not be 
confused with disc populations of BHBs. To remove BHB 
stars in the Sagittarius dwarf tidal stream from the sam- 
ple, we also eliminated stars with S < 5° . There were 298 
stars with 18.65 < 30 < 19.15, -0.3 < (g - r)o < 0.0, 
0.8 < {u — g)o < 1.5, and that passed the photometric low 
surface gravity cut. Of these, 260 had measured metallicities 
that are shown in Figure 4. The mean of the distribution is 
[Fe/H]= -1.69 ± 0.04, and the width is a = 0.71. Note that 
the error quoted here is a statistical error. Because very few 
of the SEGUE calibration stars are as blue as A stars, the 
systematic errors at this end are about ±0.3. The compari- 
son stars are about 45 kpc from the Sun, and above b = 45° , 
so they come from the distant halo. The blue straggler con- 
taminants are 10 to 28 kpc away, which is also fairly far from 
the disc at high Galactic latitudes. 
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Figure 4. Metallicity distribution of the moving group. We se- 
lected all stars with spectra from SDSS DR7.1 which had the 
same photometric properties as the moving group BHB stars 
(18.65 < go < 19.15, -0.3 < 3 - ro < 0.0, 0.8 < {u - g)o < 1.5, 
and split photometrically by likelihood of low surface gravity), 
with Galactic latitude b > 45°. The 298 stars selected are all 
likely spheroid BHB stars, of which 260 had measured metallic- 
ities. The distribution of metallicities of the spheroid BHB stars 
with similar photometric properties (thin lines) is shown along 
with the metallicities of the seven BHB stars in the moving group 
(thick lines). The moving group has a mean metallicity of —2.5, 
while the spheroid BHBs have a mean metallicity of —1.69. 



A comparison of the stars in the moving group with the 
similarly selected stars in the spheroid with a t-test gives a 
probability of 1 in 10^ that the two groups of stars were se- 
lected from the same stellar population. The Mann- Whitney 
test (also known as the Wilcoxon rank sum test) gives a p- 
value of 0.00023, or a 2 in 10,000 probability that the two 
samples are drawn from the same population. This leads us 
to conclude that the moving group is real. 

As a further test, we selected the other 8 stars in the 
lower panel of Figure 1 that have the same Vgsr and ap- 
parent magnitude as the blue horizontal branch stars in the 
moving group, and find that they have a mean metallicity 
of -1.87 with a sigma of 0.86. The t-test (p=0.948) shows 
that this distribution is not distinguishable from the back- 
ground population. (Performing the Mann- Whitney test for 
this sample is problematic since several of the stars are not 
part of our selected background distribution.) 



7 SYSTEMATICS OF DR7 METALLICITIES 
OF HALO BHB STARS 

In this section we explore the accuracies of the DR7 metal- 
licity measurements of BHB stars as a function of S /N and 
in comparison with published metallicities for globular clus- 
ters. We also explore the metallicity of BHB stars in the halo 
as a function of apparent magnitude. 

Since all of our moving group stars have low S/N, we 
need to test the accuracy of the metallicity measurements. 
We do this by selecting BHB stars in the leading tidal tail 
of the Sgr dwarf galaxy that are about the same distance 
from the Sun as the moving group. We selected all of the 
spectra from SDSS DR7 that had elodierverr> (if the ra- 
dial velocity error when matching to the stellar templates is 



10 < S/N < 25 (N=28) 
7 < S/N < 10 (N=52) 
< S/N < 7 (IM.14) 
Sgr North BHBs near gQ-19 



nn 



Figure 5. Measured metallicities of BHB star in the Sgr dwarf 
leading tidal tail, as a function of S/N. The metallicities of 
stars with S/N> 10 are considered fairly reliable. The stars with 
7 <S/N< 10 have a similar distribution. Stars with <S/N< 7 
have a significantly higher average metallicity. Therefore, we con- 
clude that the metallicity measurements of BHB stars with S/N< 
7 are not reliable. Note that the excess of stars at [Fe/HJwBG = 
—3.0 is due to edge effectes in the WBG estimation procedure, 
and does not represent a cluster of stars with identical metallicity. 



less than zero, that indicates that the object is not a star), 
0.8 < {u-g)o < 1.5, -0.3 < (ff-r)o < 0.0, 320° < I < 350°, 
30° < b < 70°, -100 < Vgsr < 100 km s'S logg9< 3.0, 
and 18.5 < go < 19.6. The 94 stars that are likely BHBs 
in the Sgr dwarf tidal stream were divided into three files 
based on S/N. There were 14 stars with 0.0 <S/N< 7, 49 
stars with 7 <S/N< 10 (excluding three with unmeasured 
metallicities), and 28 with 10 <S/N< 25. The metallicity 
distributions of these three S/N groups is shown in Figure 
5. 

The highest S/N group (10 <S/N< 25) has a mean 
[Fe/H]wBG ~ —1.94, which is similar to the group with 
marginal S/N (7 <S/N< 10) which has a mean [Fe/HwBG = 
—2.13. A Mann-Whittney test comparing these two samples 
gives a result of U — 818, z — —1.39, and a p-value of 0.165. 
These distributions are not significantly different from each 
other. This gives us some confidence that the metallicities 
of the moving group stars, which have S/N measurements 
in this marginal range, are usefully measured. 

In contrast, the lowest S/N group (0 <S/N< 7) has a 
mean of [Fe/HJwec ~ —1.42. A Mann-Whittney test com- 
paring this sample with the highest S/N sample produces 
U = 275.5, z = -2.11, and a p-value of 0.0349. This is a 
statistically significant difference, and metallicities of BHB 
stars measured by this technique are probably not accurate. 

To understand the likely systematics in our metallic- 
ity determinations, we selected BHB stars from six globu- 
lar clusters that had spectroscopic measurements of BHB 
stars in SDSS DR7. We selected stars that have the colors 
of BHB stars, are within half a degree of the centre of the 
globular cluster, and have radial velocities within 20 km s~^ 
of the published value for the cluster. M2, M3, and M13 
have metallicities very close to —1.6; M53 has metallicity of 
-1.99; and the two cluster s M92 and N GC 5053 have metal- 
licities very close to — 2.3 l|Harrisll 19961 ) ■ Because there were 
very few stars with spectroscopy in each globular cluster. 
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Figure 6. Metallicity measures of known clusters. We show the 
metalli city distribution of SPS S DR7 BHB stars, as measured 
by the IWilhelm. Beers, fc Gravl l ll999l) technique, for BHB stars 
selected from clusters of known metallicity. Uppermost panel: 
Stars in M2 ([Fe/H] = -1.62), M3 ([Fe/H] = -1.57), and M13 
([Fe/H] = -1.54). Second panel: Stars in M53 ([Fe/H] = -1.99). 
Third Panel: Stars in M92 ([Fe/H] = -2.28) and NGC 5053 
([Fe/H] = —2.29). Lowermost Panel: The seven Hercules stream 
BHBs. Note that the measured metallicities are compressed to- 
ward [Fe/H]=— 1.9 (higher metallicity stars are measured system- 
atically low, and lower metallicity stars are measured systemat- 
ically high). Also note that the Hercules stream stars appear to 
be lower metallicity than M92 and NGC 5053. 



we combined the data from clusters with similar metallici- 
ties before plotting the histograms in Figure 6. All of these 
globular clusters have horizontal branches near go — 16 (the 
magnitude range for which we have the most complete sam- 
ple of BHB stars), and thus the spectra have high S/N. 

Figure 6 shows that the SDSS DR7 WBG metalhci- 
ties for BHB stars tend to push the metallicities toward 
[Fe/H]=— 1.9. The measured metallicities are correlated 
with the published metallicities of the globular clusters; how- 
ever, stars that have higher metallicities are measured with 
metallicities that are too low, and stars with lower metallic- 
ities are measured systematically too high. 

We show for comparison the measured metallicities of 
the seven stream stars. The measured metallicity is lower 
than that of M92 and NGC 5053. 

Because the validity of our stream detection depends 
on our ability to separate the stream from the stellar halo 
in metallicity, we need to show that the metallicities in our 
background population are accurate. We selected all of the 
stars with spectra in DR7 that have photometry consistent 
with being a BHB star, as explained in §2, b > 45°, and 
5 > 5° . We further restricted the sample by insisting that 
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Figure 7. Metallicity of the outer halo vs. distance. We show 
the metallicity distribution of SDSS DR7 BHB stars in several 
apparent magnitude ranges. The last panel shows the six Hercules 
stream BHBs (thick lines in lower left of diagram) that survive the 
additional cut in logg. They arc inconsistent with being drawn 
from the same distribution as the other stars in this panel. All 
of the stars have ugr colors consistent with BHB stars, b > 45° , 
5 > 5°, and logg < 3.75. The [Fe/H] and logg measurements 
are from the WBG llWilhelm. Beers, fc Grav|ll999h techniques, as 
implemented in the SSPP of DR7. A Gaussian with mean —1.9 
and sigma 0.45, normalized to the number of stars in each panel, 
is shown for reference. The mean metallicity does not change as 
a function of distance from the Sun. The distribution is wider 
in the lowest panel because the signal-to-noise is lower here. If 
one selects only the highest S/N spectra in this bin, the width is 
similar to the four other panels. 

the WBG estimate of \ogg, as measured in DR7, was less 
than 3.75. 

In Figure 7, we show the DR7 WBG metallicity as a 
function of apparent magnitude. We have spectra for stars 
with 14.5 < go < 19.15, which span distances of 6 to 50 kpc 
from the Sun, all at high Galactic latitude. The metallicity 
distributions for all of the stars brighter than go < 18.5 
are similar, with a mean near [Fe/H]=— 1.9 and a sigma 
of 0.45. In the faintest set of stars, which have apparent 
magnitudes similar to that of the newly detected stream, 
the distribution appears considerably broader, but with a 
similar mean. We attribute the increased width to the lower 
S /N in many of these spectra. The mean is somewhat lower 
in the last panel in Figure 7 than in Figure 4, probably due 
to a cleaner sample of BHB stars, with less contamination 
from BS, from the additional restriction in loggr. 

We show for comparison the six stream stars that have 
logg < 3.75 in the last panel of Figure 7. The metallicity 
distribution of these stars is inconsistent with the distribu- 
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Figure 8. Projection of the moving group in equatorial coordi- 
nates, as detected in SDSS DR7. The small black dots show the 
positions of A stars with spectra that are not from SEGUE plates 
as selected from DR7. Only twelve of these stars (large cyan sym- 
bols) have surface gravities, metallicities, apparent magnitudes, 
and line-of-sight velocities that are consistent with the moving 
group, as detected in Figure 9. Of these twelve, four are very 
tightly clustered and a fifth is only a few degrees away and in line 
with the other four. 



tion of the other stars in this figure at the 98% confidence 
level, as determined by the Mann- Whitney test. 

The value that we find for the metallicity of the outer 
halo, as measured from SDSS BHB stars is on the low side 
of the di stribution of p reviou s measurements, but higher 
than the ICaroUo et alj |20o3) metallicity measurement of 
[Fe/H]=-2.2. (We note, however that CaroUo et al. 2007 
distribution of BHB stars in their paper, and the BHB star 
distributions peak closer to [Fe/H]=— 2.0.) 




Figure 9. Detection of the moving group in SDSS DR7. This 
figure is similar to the upper right and lower left panels of Figure 
1, with the additional restriction that the all the stars in these 
panels have [Fe/H]wBG < —1-95. The significance of the group 
over background is enhanced with this metallicity cut. 



8 THE MOVING GROUP IN DR7 

The moving group was originally discovered by searching 
through data from SDSS DR6. In this section we show the 
moving group as selected from the now public SDSS DR7 
databases, and using our knowledge of metallicity and sur- 
face gravity determination at low S/N, as learned in the 
production of this paper. 

We selected from SDSS DR7 all of the spectra with 
elodierverr > (which selects stellar spectra), magnitude 
15 < po < 23, the colors of an A star [0.8 < {u — g)o < 
1.5,-0.3 < ((? — r)o < 0.0], and a flag that indicates the 
spectrum is not from a SEGUE plate. We did not select 
stars from SEGUE plates since most of them have not been 
searched for moving groups of A stars, and because their 
non-uniform sky coverage, different section criterion, and 
different exposure lengths make it difficult to use them in 
a statistical calculation. We show the distribution of SDSS 
DR7 A stars in equatorial coordinates in Figure 8. 

We then selected the low surface gravity (loggO < 3.0) 
stars with 50° < I < 80° and 40° < 6 < 55°, to match 
the region of the sky presented in Figure 1, and metallicity 
[Fe/H]wBG < —1.95, since we know we are looking for a low 
metallicity moving group. Figure 9 shows plots like those 
used in the moving group discovery, but using data from 
this low metallicity DR7 dataset. From the upper panel, we 



identify a clump of stars with —35 < Vgsr < 0.0 km s"'^ and 
18.8 < go < 18.9. The lower panel shows the positions in 
Galactic coordinates of the same objects from the top panel, 
with filled circles showing the six objects that are clumped in 
line-of-sight velocity and apparent magnitude. Four of these 
objects are very tightly aligned in angular position on the 
sky. 

To get a sense for how common a moving group with 
these characteristics is, we selected from our original DR7 
dataset all of the stars with surface gravity, metallicity, line- 
of-sight velocity, and apparent magnitude that are similar to 
the stars in the clump. In the entire northern sky portion of 
the SDSS, there are twelve stars with these characteristics. 
Four of the stars are in a tight line on the sky, and a fifth is 
only a few degrees away from the clump along the same line. 
Because this fifth star is a candidate member of the moving 
group, and not included in the original list of seven, we have 
added it to the bottom of the list of candidate members in 
Table 1. The four highly probable members of the moving 
group, which were selected from both the DR6 and DR7 
datasets, are marked in Table 1 with a footnote. 

The four highly likely moving group members are a sub- 
set of the original seven, and have velocity and metallicity 
characteristics that are nearly identical to our previous re- 
sults. The average Vgsr is — 10±3 km s~^, with a sigma of 6 
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km s~^. The average metallicity is [Fe/H]wBG ~ —2.4 ±0.2 
with a sigma of 0.39. 

Following a similar procedure to that used in Section 
5, we now calculate the probability that this group of four 
stars is a random coincidence. There are 6118 stars with 
b > 35° in the A star sample depicted in Figure 8. Of these 
stars, 3335 have -9 < logg9 < 3.0 (54.5%), 1171 have -9 < 
[Fe/H]wBG < -1.95 (28.9%), 224 have 18.8 < go < 18.9 
(3.7%), and 487 have -23 < Vgsr < km s"! (8.1%). The 
four stars in the clump have sky positions within 63.72° < 
/ < 66.62° (0.8% of longitude range) and 47.37° < 6 < 
48.92° (4.2% of area-corrected latitude range). Because the 
density of the stars in Figure 8 is not uniform, we need to use 
a multiplier on the fraction of stars in the longitude range, 
since that is fractionally the smallest dimension. In the 3.0 
square degree area that contains the four stars, there are 
8 stars in the original sample of 6118. So, in the region of 
the clump there are 2.7 stars per square degree, while there 
are 6118/8796 = 0.69 stars per square degree in the whole 
sky region, a factor of 3.9. Using the procedure outlined in 
Appendix A, we calculate that we would expect 0.0038 stars 
in a 6-D region that instead has four stars. The expected 
number of clumps of four stars is 0.0238. The p-value is 
therefore less than 0.024 (p < E = 0.238). Since p < 0.05, 
this is a significant detection. 

Note that this calculation is extremely conservative. 
First, the calculation assumed that we searched all 6118 
stars for clumps, when in fact we only searched the ones that 
had low surface gravity, as determined by the photometric 
indices. Second, the four stars are not randomly distributed 
within the {I, b) box, they are tightly constrained to a line, as 
one would expect for a tidally disrupted globular cluster; if 
our statistical analysis allowed a non-axis aligned region, the 
computed statistical significance would have been stronger. 
And third, when we calculated the density of stars in the 
region containing the moving group we included the moving 
group stars (half of the sample of eight) in the density cal- 
culation. This artificially inflates our Galactic latitude mul- 
tiplier. To estimate the effect of just this third factor, we 
counted the number of A stars in a slightly larger region of 
the sky with 60° < I < 66.63° and 47.37° < b < 48.92° 
and found 16 stars, of which 4 are members of the moving 
group. This gives us a local density of 12/6.87 — 1.75 stars 
per square degree, which is a factor of 2.5 times the average 
star density. Using this multiplier, we expect to flnd 0.0025 
objects where we actually find four, and the expected frac- 
tion of the time a clump of four stars in this size region of 
6-D space should be found in this size data set is not more 
than 0.0063. 



9 ESTIMATED LUMINOSITY OF THE 
PROGENITOR 

We estimate the properties of a progenitor of this moving 
group, making the assumption that it was a globular cluster, 
and that we have detected the majority of the BHB stars 
from what was the core of the star cluster. We expect the 
progenitor was a globular cluster because the stars are very 
well aligned in the sky, have a velocity dispersion consistent 
with instrumental errors, and there are very few stars iden- 
tified. If the progenitor was a dwarf galax;y, so that it also 



has younger or more metal rich components, or if we have 
found only a knot of increased stellar density along a longer 
tidal debris stream, then the inferred size of the progenitor 
is larger. 

We detected 4-8 BHB stars in the moving group. In the 
color-magnitude range in the vicinity of the moving group, 
we have spectra of 50 of the 73 candidate BHBs (68%). 
Therefore, we expect that there a re ~ 10 BHBs in the mov- 
ing group. iBorissova et al.l ||1997D found five BHB stars in 
the globular cluster Pal 13 ([Fe/H]=-1.7 4, My = -3.74, 
Harris 1996). lYannv. Newberg et all (|2000h find at least five 
BHB stars in the globular cluster Pal 5 ([Fe/H]=-1.41, 
Mv = —5.17, Harris 1996). The moving group is consistent 
with a progenitor that is like one of the smaller globular clus- 
ters of the Milky Way, with a total integrated luminosity of 
about Mv ~ -4. 

To understand why this moving group could not be 
identified from photometry alone, we estimated the num- 
ber of red giant and subgiant branch stars in the moving 
group by comparison with NGC 2419. Note that at 50 kpc, 
turnoff stars are at go ~ 22.7, which is close to the lim- 
iting magnitude of the SDSS, so deeper photometry is re- 
quired to clearly detect significant density enhancements ex- 
pected from turnoff and main sequence stars fainter than 
these limits. We selected SDSS DR6 stars within about nine 
arcminutes of the centre of NGC 2419 (this does not in- 
clude stars near the very centre of the GC, since individ- 
ual stars are not resolved there). BHB stars were selected 
with 20.2 < So < 20.6 and -0.4 < {g - r)o < 0.0. Gi- 
ant stars were selected within a parallelogram with ver- 
tices [{g-r)o,go] = [0.5,20], [0.9, 18], [0.7,20], [1.0, 18]. Sub- 
giant stars were selected within a triangular area with ver- 
tices [ig-r)o,go] = [0.5,21], [0.1,23], [0.6,23]. We found 101 
BHBs, 98 red giants, and 302 subgiants over background in 
the region of sky that was searched. Therefore, we expect 
about the same number of red giant stars as BHBs in our 
moving group, and about three times as many subgiants. We 
shifted the magnitudes of the color-magnitude boxes by 1.41 
magnitudes (the difference in distance modulus between our 
moving group and NGC 2419) and counted the number of 
stars in 60° < / < 70° and 45° < & < 50° . The moving group 
is expected to have 10 of 73 candidate BHB stars (1 sigma 
fiuctuation), 10 of 2742 candidate giant stars (0.2 sigma), 
and 30 of 13,182 candidate subgiant stars (0.3 sigma). The 
background star counts are much too high for us to find this 
moving group in photometry alone. 



10 CONCLUSIONS 

We have detected a moving group of at least four BHB 
stars in the corner ("toe") of the Hercules constellation 
spilling into Corona Borealis. These stars are coincident in 
angular position [(/,&) — (65°, 48°)], apparent magnitude 
{go = 18.9), line-of-sight velocity {Vgsr = —10 km s~^, 
(yvggj. < 10 km s"'^, (Vr) = —157 km s~^), and metallic- 
ity ([Fe/H]=— 2.4). We expect that the progenitor of this 
moving group was a low metallicity globular cluster, with a 
luminosity like that of one of the smaller globular clusters 
in the Milky Way halo. 

We show that useful surface gravities and metallicities 
are measured for BHB stars with S/N> 7 in SDSS DR7. The 
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mean metallicity of BHB stars in the outer halo is similar to 
M53, which has a published metallicity of [Fe/H]=-2.0. The 
metallicity does not appear to change with distance from the 
Sun (6 < i? < 55 kpc). Our measurement o f the spheroid 
metallicity is slightly higher than claimed bv lCarollo et alj 
l|2007l ) and somewhat lower than earlier studies of outer halo 
stars. 

The Hercules moving group is one of many tidally dis- 
rupted stellar associations expected to comprise the spheroid 
of the Milky Way and could not have been identified from 
photometry alone; more complete spectroscopic surveys are 
required to identify the component spheroid moving groups, 
and determine the merger history of our galaxy. We present 
a statistical technique that allows us to estimate the signif- 
icance of clumps discovered in multidimensional data. 
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APPENDIX A: THE STATISTICS OF CLUMPY 
DATA 

Al Introduction 

For this paper, we needed to estimate the probability of find- 
ing seven BHB stars in a restricted portion of l,b,Vgsr, go 
parameter space. In this appendix, we present a solution 
to the more general problem of determining the probabil- 
ity of finding K oi N data points in a small portion of a 
d-dimensional space. 

Suppose one analyses of a set of A'^ data points in 
a large d-dimensional hypercube x ... x 

and finds a small, axis-aligned, d-dimensional hypercube 
[ai,6i] X ... X [ad,bd\ that has K ^ 2 data points in it. 
The natural question is whether this "clump" of K data 
points is statistically significant, or whether even random 
data would yield such a clump. Here we precisely esti- 
mate an i5-value for the expected number of clumps, thus 
providing an upper bound to the probability that random 
data would somewhere have K or more data points within 
an axis-aligned, d-dimensional hypercube with dimensions 
(wi, . . . , Wd) = {bi-ai,...,bd- ad). 



A2 Theory 

We start with the simple case where a full data set of TV 
points lies within a d-dimensional unit hypercube [0, 1]**. 
We further assume that a random (or null) model would 
distribute these N points uniformly within the hypercube. 
Additionally, wc assume that the hypercube is not peri- 
odic/toroidal. We relax these simplifying assumptions in 
subsequent sections. 

We say that an axis-aligned, d-dimensional hypercube 
contained within [0, 1]"* is a box . We say that a set of K ^ 
2 of the A*' points is a K-boxed set if there exists a box 
that contains the K points without containing any of the 
remaining N — K points. We define the minimal box for a 
/sT-boxed set as the intersection of all boxes that contain the 
K points. 

A box of dimensions (toi, . . . , Wd) or smaller with K or 
more points exists only if there is a i^-boxed set that hcis a 
minimal box with dimensions not greater than (wi, . . . , Wd). 
Thus, we aim to compute the expected number of the latter. 
We use as an integrand the probability density that a box 
will be a minimal box for a set of K points, and we integrate 
over all applicable boxes. 

Consider a box x ... x [/d,3d] with dimensions 

[li, . . . ,ld) = {h - gi, ■ ■ ■ , fd - Qd) and hypervolume 



(Al) 



It will be a (non-minimal) box indicating that a set of K 
points is a if -boxed set if exactly K points fall within it. 
Thus the probability that it indicates a if -boxed set is 



(A2) 



However, it is our goal to calculate the probability density 
for a minimal box for K points. This box is minimal only if, 
for every dimension index i, we have that /, is the minimum 
of the K points' ith coordinate values and Qi is the maximum 
of these coordinate values. Thus, for each dimension, we 
must have one of the K points achieve the minimum and 
one of the remaining K — 1 points achieve the maximum. 
The remaining K — 2 coordinate values can be anywhere 
in the range {fi,gi)- (Note that the possibility of having 
multiple points exactly achieve the minimum or maximum 
is an event that has a probability density of zero and thus is 
safely ignored.) Thus, the ith dimension contributes a factor 
of K{K— to the probability density. The probability 

density for the event that the box is a minimal box for a K- 
boxed set is 



, Id 



K'^{K-l)'^V'^-''{l-Vy 



.(A3) 



We note that a box with dimensions (Zi, . . . , Id) can have 
its minimum corner {fi, ■ ■ ■ , fd) anywhere in [0,1 — h] x 
. . . x [0, 1 — Zd]. When we restrict attention to minimal boxes 
with dimensions not greater than {wi, . . . , Wd) the expected 
number of if -boxed sets is 



E{K, wi, ...,Wd) 



/ dh... did[Y[ 

Jo Jo 



{1-h) Pr[if|«i,...,Zd](A4) 
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Equation IA4I provides an upper bound to the proba- 
bility that random data will exhibit a clump of K or more 
points in a box of dimensions (wi, . . . , Wd) or smaller, be- 
cause such a box exists only if there is a if-boxed set 
that has a minimal box with dimensions not greater than 
{wi, . . . , Wd)- Under many circumstances, when E <^ 1 the 
probability of multiple A'-boxed sets is small and the bound 
is quite tight. 

A3 Implementation 

We observe that 

N-K 



N-K 



(A5) 



and, to speed the integration using Maple 7, we seek a way 
to truncate the sum at significantly fewer terms. We observe 

that ^ y™ is roughly a bell curve as a function of 

m. The bell curve has mean and standard deviation given 
by 



{N -K)V 




(A6) 



(A7) 



where we have assumed that K ^ NV , i.e., that a clump 
is under consideration only when the number of points it 
contains is more than the expected number. 

Thus, we can include all terms of (1 — 1/)^"^ out to 
z = IQ standard deviations with m* = max([/^ -I- zct],25) 
and the precise approximation 



(A8) 



We use this approximation in Equation I A4I and we evaluate 
the integral with Maple 7. 

A4 Variations 

A^.l Non-Unit Hypercubes and Non-Uniform Random 
Distributions 

In many circumstances, including the case presented in this 
paper, the set of A'^ data points does not fall uniformly into 
the unit hypercube under the random/null model. For in- 
stance, one does not expect the radial velocities of spheroid 
stars to be uniformly distributed in any velocity range, and 
they are certainly not limited to the range [0, 1]. 

We say that the density is separable if it can be writ- 
ten as the product of functions, each of which depends on 
a single dimension Pr(x) oc Y[i^i^i{^i)- For instance, the 
uniform spherical density, 

Pr(r, e, <P) = -^^^^r^ sm{e) dr dd d^ . (A9) 

is separable. When a probability density for a coordinate 
vector X G [yli,_Bi] x ... x is separable, we can 

perform a coordinate transformation to a coordinate vec- 
tor y that is uniformly distributed in the unit hypercube 
and which preserves axis-aligned hypercubes: 



I a' hiixOdx^ 



(AlO) 



The clumping question then becomes one of percentiles, i.e., 
what is the probability that random data will show K or 
more points in a axis-aligned hypercube with dimensions 
that do not exceed Wi percentiles x . . . x w^j percentiles. In 
this paper, the radial velocities of the stars and the apparent 
go magnitudes were treated as separable dimensions. 

A4.2 Periodic or Toroidal Boundary Conditions 

In many circumstances, it is desirable to allow Jf-boxed sets 
to wrap around the boundaries in one or more dimensions. 
In this paper, it was desirable to allow a box to span the 
prime meridian in a longitude coordinate. In such situations 
the restriction that the minimum corner (/i, . . . , fd) fall in 
[0, 1 — /i] X . . . X [0, l~ld] is relaxed. When the ith coordinate 
is periodic, fi can be any value in [0, 1]. Equation I A4I is thus 
modified by removal of the (1 — h) factor for each periodic 
dimension i. 



A4.3 Poisson Distributions 

In many circumstances, it is desirable to have the ran- 
dom/null model be Poisson; instead of requiring exactly A'' 
points in the unit hypercube, the number of points n is se- 
lected using a Poisson distribution with mean A'^: 



Pr[7il7V] 



(All) 



In our example, the Poisson distribution applies if there are 
an average of 4149 BHB stars with spectra in the region 
6 > 45°, whereas the uniformly random model would put 
exactly 4149 BHB stars with spectra in that same region. 
In the Poisson case. Expression IA2I becomes 

(A12) 



N 



V e 



Equation IA3I becomes 

PT[K\h,...,ld] 



N'' ^ 



(A13) 



and Equation I A4I is unmodified, except that it incorporates 
Equation lA 131 instead of Equation I A3I 

For integration, instead of truncating the series for (1 — 
we truncate the series for 



E 

m=0 



(A14) 



We observe that (NV)^ jm}. is roughly a bell curve as a func- 
tion of m. The bell curve has mean and standard deviation 
given by 



A* 



(A15) 
(A16) 



Thus, we can include all terms of exp(— A"!/) out to z stan- 
dard deviations with m* = ma.x{\ ij,-\-za~\ , 25) and the precise 
approximation 



E 

m = 



i-Nvy 



(A17) 
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We use this approximation in Eg nation I A 1 31 and we evaluate 
the integral of Equation IA4I with Maple 7. We determined 
that for this paper the distinction between random and Pois- 
son distributions was unimportant. 



A4-4 N on- Separability 

When the dimensions are not separable, the assumption that 
the joint density function is the product of the individual di- 
mension's density functions does not hold. However, if there 
is a set of individual density functions for which the assump- 
tion is approximately true, we define the overdensity to be 
the ratio of joint density function to the value implied by 
multiplying the approximating density functions together. 
The calculations for the separable case make a good approx- 
imation for the non-separable case if the percentile widths 
{wi, . . . ,Wd) are adjusted by factors whose product is the 
overdensity. When the overdensity is greater than 1.0, a con- 
servative approach is to multiply the smallest of the widths 
of the non-separable dimensions by the overdensity. 

This paper has been typeset from a Tff)i/ I^Tp^jX file prepared 
by the author. 



